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Abstract 

The scaling function for the critical specific heat is obtained exactly 
for temperatures above the bulk transition temperature by working in 
the spherical limit. Generalization of the function to arbitrary a (the 
specific heat exponent), gives an excellent account of the experimental 
data of Mehta and Gasparini near the superfluid transition. 
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One of the primary theoretical tools for a strong coupling problem is the 
spherical limit and the ^ expansion [1-5]. In critical phenomena, it has often 
provided useful insight in the past and is still one of the first techniques to be 
tried when a new kind of critical phenomenon surfaces [6-11]. A problem that 
has been receiving a fair amount of experimental and theoretical attention is 
that of the specific heat of liquid helium near the lambda point in a confined 
geometry. This is a very effective way of studying the finite size effects 
(FSE) since by going close enough to the lambda point (particularly easy in 
a space-shuttle experiment), the correlation length can be made larger than 
the confining length if the liquid is enclosed in a thin wafer geometry. The 
bulk specific heat C p diverges near the lambda point as ^ a l v , where £ is the 
correlation length proportional to | T — T c \~ u . For the finite size system, 
the specific heat remains finite at T = T c , with the value proportional to 
L a / U where L is the confining length. For a given L, the specific heat changes 
from a £ a / u to a L a l v as one approaches the bulk lambda point. The change 
from one regime to the other is described in terms of a scaling function f(x) 
in terms of which we can write the specific heat as 

C p (£, L) = C^f{i/L) + Constant (1) 

where the function f(x) has the property that /(0) = 1 and f(x) ~ x~ a / u 
for x 1 . It is the function f(x) that has been the object of several 
experimental and numerical investigations. In this letter, we determine f(x) 
in the spherical limit. We also show how the result can be put to practical 
use. 

We consider a n-component Ginzburg-Landau model with the free energy 
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functional 

F = / ^ E <^ E (V0,) 2 + A (E ] (2) 

i=l i=l i=l 

where m 2 ~ (T — T ) , T being a transition temperature and the coupling 
constant A is 0(n~ l ) so that the quartic term remains the same order as 
the quadratic term when n — > oo (the spherical limit). We consider the 
system confined in one direction (we will call this the z - direction) within 
an extension L and the boundary conditions will be taken to be of Dirichlet 
type , i.e., fa = at z = and z = L . The expansion in Fourier modes for 
fa(r) is 

/ — * \ ^ f d k ^ — ^ . — * — * Tin ' z . . 

<M r ^) = "7= / /„ xni ^fan\K)VXVlk.Rsin—- (3) 

The specific heat is given by the correlation function 

C p = JlJ d D - 1 r 1 d D ~ 1 r 2 dz x dz 2 < E $(rl, z x ) E 4>)(r 2 , z 2 ) > 

i j 
1 r d D - 1 p 1 d D - 1 p 2 ^ ^ 

= 4L / (2^ < ? 

E fa™ (Pi) 0jn 2 (-P 2 ) > (4) 
j 

The free energy of Eq(2) written in terms of the fields fa n because 

p = i:S£^\\^ 2 +k 2 +^)fa n {k)fa n {-k)} 

A vv r d D ~ l h d D ~ l k 2 d D ~ l h , 
+ A Z^2^J ( 2vr )fl-i (2n)D-i (2tt) d - 1 

faih, z) fa(k 2 , z) fa(k 3 , z) fai-h -k 2 - k 3 , z) (5) 

and the averaging shown in Eq(4) has to be done with this free energy func- 
tional. 
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The Gaussian limit has to be disposed of first i.e. the limit where A = 
and the action is quadratic. 

< 0i,m(Pl) <Pj,n 2 {P2) > 



AT /■ 1 



2L ^7 (27r) D - 1 (p2 + m 2 + 2 ^) 

NL 4 /■ rf D -V r coth^p 2 L 2 + m 2 L 2 



7 f27r) D - 1 [ 



2L 7 (27T) 15 - 1 L 2(p 2 L 2 + m 2 L 2 ) 3 / 2 

I cosech 2 y/p 2 L 2 + m' 2 L' 2 



(p 2 L 2 + m 2 L 2 ) 2 2(p 2 L 2 + m 2 L 2 ) 

— Jjj(mL) (6) 



where 



I D {mL) 



d p coth^/p 2 L 2 + m 2 L 2 



(2tt) d - 1 l 2(p 2 L 2 + m 2 L 2 ) 3 / 2 

1 cosech 2 Jv 2 L 2 + m 2 L 2 , 

+ o^ 2 ; 2 , ] (7) 



(p 2 L 2 + m 2 L 2 ) 2 2(p 2 L 2 + m 2 L 2 ) 

For L 3> this gives the usual gaussian limit answer C p ~ m~( 4 ~ D ^ 2 , 
while for m — > at finite L, C p ~ £( 4 --°)/ 2 . The scaling function is provided 
by J D . 

We now turn to the spherical limit. The propagator < <pi n (k) 4>i n (—k) > 
has the structure (M 2 + A; 2 + ^^r-) -1 where M is the renormalized mass 
dressed by the bubble [4] in the spherical limit. As usual M ~ (T — T c )~ v ', 
where in the spherical limit v = (D — 2)~ 1 . The specific heat graphs [4] for 
the spherical limit are shown in Fig.l. The first graph (fa) corresponds to 
the gaussian limit. Each additional loop brings an interaction of strength 
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iV -1 and a combinatoric factor N (for large N). Thus, each of the graphs 
(lb, lc, Id. ...etc) are the same order as the single loop and the sum defines 
the spherical limit. For L — > oo, the contributions factor and the successive 
contribution factors are I 2 (lb) ,/ 3 (lc)....etc. The total specific heat is 
proportional to y^j . In the case of finite geometry this simple geometric 
series does not obtain as we show below. 

To understand the complication, let us look at the two loop graph. The 
contribution from the graph is , 

= Ej^^^ik < ^(Pl) hnA-Pl) f d* E 

E 4>s, mi (ki) 4>s, m2 (h) 4>t, ms (h) 4>t, m4 (-h -k 2 - h) 

mi ,m 2 ,m3 ,1714 

miiTZ . m 2 nz . m 3 nz . m^nz . . 
sin—^— sin—^— sm—j— sin—j—<p jin2 {p 2 ) <Pj,n 2 {-P2) > 

, AT 2 /■ ^-Vi ^ D ^V 2 r , / \ 1 / x i2 

= U J (2tt) d ~ 1 (2tt) d - 1 [ < ^ n ^ Pl) ^(-Pi) > \ 
[ < (f>j,n 2 {p 2 )(f>j,m 2 {-P2) > J y — — sm — — 

= _ 4^nS^ [ < ^' ni(Pl) 0i ' ni( ^ l} > ] 

r / / \ / / \ 12 A 7 /■, 2niiiz 2n 2 nz 
[< <Pj,n 2 {P2) <Pj,n 2 {-P2) >\ J d,Z (1 + COS — ) (1 + COS- 



N r d^ Pl d^p 2 ^ 



L 

m,n 2 ' 

< <Pi,n 2 (p 2 ) <f>i,n 2 (-P2) > ] 2 ( 1 + ^6 r 



V2f 1 1 • 

1 ~ L 2" n l' ra 2 



O L 7 D + , r2 2^ 



2 1 " 4L 2 V J (I*) - 1 (2tt) d - 1 

1 



(p2 + m 2 + ^) 2 (p 2 + m 2 + ) 2 j 

+ (8) 



where 

/<i> = J_V f d °~ lpi d °' lp2 1 1 (q) 

D 4L*^J (27T)*-1 (2n)D~l ( p 2 + M 2 + )2 (p 2 + M 2 + n^ )2 ^ J 

The difficulty is obvious from Eqn(8). The presence of the term J$ means 
that the series for the specific heat is not a geometric series. The three loop 
term analysed in a similar manner yields f {II + 2I D J$ ) + , where 



' D (2L) 2m ? / (2tt) d - 1 (2n) D ~ 1 "' (27T) - 1 

1 1 



•J , 



(p\ + + S^L)2 (pi + MH ^ (pl +l + + 

(10) 

The four loop calculation is -f [i£ + S/^J^ + 2/ D jg } + j£ )2 + }. 
Similarlythe five loop term is f [I^+AIfjJ^ + 3I D J^ )2 +2I D J { ^ +2J { ^ jff + 
J^J.The contribution from the six loop term is — y [if, + hl^J^ 1 + 2Ipjff + 
6/|, J« 2 + 6/ D J« J( 2) + J« 3 + 3/£jg> + 2/ D jg> + 2j£>jg> + Jg )2 + 4 5) ] 
and so on. 

The pattern is now clear. The specific heat is obtained as 
N 



C = -j{[Id-I 2 d + I S d + }-[J^-J^ + J^-.][l-2I D + 3ll-..} 



r (l) j(2) j(3) 

~ J% ] + 4? " -.] 2 [1 " 3/z, + 6/ 2 - ..] - [J« - J« 
+J§ ) -..] 3 [1-4/d + lOTf, -..]....} 

r -^D Em^fl ^ (Em (Em^j) _ i 

2" 1 lT7^" (1 + /d) 2 (1 + /d) 3 (1 + /d) 4 " j 

AT T ^ fH ^ tH ^ T(m) 

2 4 + (1 + /d) 2[ (1 + / d ) + 1 



r(m) 

J ' r J K L~im J D ^ \ 

~2 i l + I D ~ (l + l D ) 2 , 1 

1 (Wo) 
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N Id — Em Jp^_ q-jn 
2 l if tM ^ ' 

This is the final answer for any dimension. The first part of Eqn(lO), where 
we generalize the pattern can be proven by induction [12] . We assume that 
this form is consistent with the form obtained after a m-loop calculation and 
prove in a long but straightforward fashion that it holds for the (m + 1) loop 
as well. We now write down the explicit expression for I D and in D = 3 

[ 'CothML-^-} (12) 



° 2M 1 ML 
For L — > oo , J3 = M" 1 as expected while for M — > , 1$ — j in accordance 
with finite size scaling approximations. Straightforward algebra yields 



and for general m 



3 oCo™\™l r2m 7l/r2m V / 



2(2m)m! L 2 ™ <9M 2 
which allows us to write Eqn(lO) as 



LF(JM 2 L 2 + \) 
1 + LF(,/M 2 L 2 + |) 



where 

F(X) = 1 [ Co^X - 1 ] (16) 

The spherical limit scaling function is contained in Eqns(14) and (15) - the 
central result of our paper. 

We now show how the above equations can be used in the practical situa- 
tion of liquid He 4 near the superfluid transition in D = 3 . The experimental 
specific heat is known to be almost logarithmic, while in the spherical limit 
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in D = 3 , ^ = — 1 . We can split off a background part from Eqn(14) by 
writing 

c = ! 1 

i + lf(^Jm 2 l 2 + 1) 

(17) 



LF(y M 2 L 2 + \) 

The second step follows from the fact that in the critical region the contribu- 
tion LFiJ M 2 L 2 + 7) is expected to be large. It is the factor , 1 

which describes the scaling function with - — —1 . If we want to write the 
answer in terms of ^ , then the specific heat is [ LF(^JX 2 + |) ] ^ and for 
- — > as is the experimental situation, we can write 



C(M,L) = C In [ LF(J M 2 L 2 + ^) ] + constant 



C In [ . AL (Coth x M 2 L 2 + 



\JM 2 L 2 + I V 4 yjM 2 L 2 + \ 

(18) 



where A = k^o- F° r L — > 00 , the bulk specific heat 

C = C ln[±] 
= Co In I 



K t v 

= vColn^] (19) 

The experimentally measured bulk value has this form [13] and we identify 
vCo = 5.3 and t = \ . The value at M = is C In [2KL(Coth\ - 2)]. 
Subtracting the C In AL from Eqn(18) to obtain a pure scaling part 

AC(M, L) = C(M, L) - Co In AL 



where X = ML . We show this plot in Fig(2) and the data of Mehta and 
Gasparini [14]. The agreement is the proof of the practical utility. 
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Figure caption 

Fig.l. Specific heat graphs in the spherical limit are shown. 
Fig.2. AC = C(M,L) - C lnAL ,Eqn(20), plotted against (ML) l ' u . The 
solid curve refers to our theory and the data is taken from the experiment of 
Mehta and Gasparini [14] 
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